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ABSTRACT 

It has now become clear that the radio jet in the giant elliptical galaxy M87 must turn on very close to the black 
hole. This implies the efficient acceleration of leptons within the jet at scales much smaller than feasible by the 
typical dissipative events usually invoked to explain jet synchrotron emission. Here we show that the stagnation 
surface, the separatrix between material that falls back into the black hole and material that is accelerated outward 
forming the jet, is a natural site of pair formation and particle acceleration. This occurs via an inverse-Compton 
pair catastrophe driven by unscreened electric fields within the charge-starved region about the stagnation surface 
and substantially amplified by a post-gap cascade. Eor typical estimates of the jet properties in M87, we find 
excellent quantitive agreement between the predicted relativistic lepton densities and those required by recent high- 
frequency radio observations of M87. This mechanism fails to adequately fill a putative jet from Sagittarius A* 
with relativistic leptons, which may explain the lack of an obvious radio jet in the Galactic center. Finally, this 
process implies a relationship between the kinetic jet power and the gamma-ray luminosity of blazars, produced 
during the post-gap cascade. 

Subject headings: accretion, accretion disks - black hole physics - galaxies: individual (M87) - galaxies: jets - 
gamma rays: galaxies - radio continuum: galaxies 


1. INTRODUCTION 

Roughly 10% of active galactic nuclei (AGN), exhibit pow¬ 
erful radio jets. These can extend for intergalactic distances, 
and potentially have impacts for the formation and evolution 
of their galactic environments. Great strides have been made in 
the past two decades in the theoretical modeling of the structure 
and formation of AGN jets, with current general relativistic 
magnetohydrodynamic (GRMHD) models now capable of re¬ 
producing the structure and Lorentz factors (L) of observed ob¬ 
jects. In these, large-scale electromagnetic fields sourced and 
collimated by the surrounding accretion flow efficiently extract 
the rotational energy of the black hole (Blandford & Znajek 
1977; Tchekhovskoy et al. 2010,2011, 2012; Tchekhovskoy & 
McKinney 2012; McKinney et al. 2012; Tchekhovskoy 2015). 
Within the context of GRMHD simulations, a canonical out¬ 
flow structure has emerged, consisting of a central relativistic 
jet core surrounded by a moderately relativistic jet sheath, fol¬ 
lowed by a non-relativistic disk wind. 

Less certain is why jets shine. This has become especially 
pressing now that it has become clear that the jet in M87 ex¬ 
hibits horizon-scale structure, implying emission on horizon 
scales (Doeleman et al. 2012). It is widely believed that from 
the radio to the UV, the observed jet emission is due to syn¬ 
chrotron emission from a population of nonthermal leptons. 
Not presently known is the mechanism by which this popula¬ 
tion is produced. Uncertain are the processes by which jets are 
mass loaded, particles within the jet are accelerated, and the dy¬ 
namical structures associated with the observed jet at different 
wavelengths (e.g., the jet spine vs. jet sheath). Any candidate 
mechanisms must simultaneously explain two critical features. 

First, the leptons must be located within the fast region of 
the jet. Radio observations far from the radio core show super¬ 
luminal motions which imply large values of the bulk Lorentz 
factor, r 1. Thus the synchrotron emitting particles in par¬ 


ticular must be moving at these extreme velocities. This sug¬ 
gests that the radio emission is not due to a mass-loaded, trans- 
relativistic disk wind surrounding and potentially collimating 
the jet. This is problematic because jet formation simulations, 
as well as approximate analytical models, imply the presence 
of large-scale, dynamically dominant, ordered magnetic fields. 
Correspondingly, the exchange of particles onto the jet mag¬ 
netic held lines is dramatically suppressed, with the conse¬ 
quence that this is unlikely to be the mechanism by which jets 
are mass loaded. Funnel regions can potentially be filled with 
pairs via collisions of gamma-rays formed in the disk. These 
gamma-rays are presumably produced through a sequence of 
inverse Compton scatterings of disk sub-mm photons with the 
hot disk electrons (Levinson & Rieger 2011; Moscibrodzka 
et al. 2011), though such a mechanism requires order unity ra¬ 
diative efficiency of the accretion flow. However, for the ra- 
diatively inefficient accretion flows believed to be relevant at 
the low Eddington ratio observed in objects like M87, most of 
the sub-mm emission is produced in the jet, and this process 
is suppressed. Alternative processes, such as neutron diffusion 
(Levinson & Eichler 2003), have been suggested as a means 
to accomplish this, though fail at low source luminosities (e.g., 
M87). 

Second, the leptons must be accelerated. Simply mass load¬ 
ing a jet is insufficient; assuming the primary emission mech¬ 
anism is synchrotron emission, the leptons must also be ac¬ 
celerated to large internal Lorentz factors as well. Far from 
the black hole this may be accomplished via internal shocks, 
magnetic reconnection, damping of plasma waves, etc. How¬ 
ever, already within lOrg bright radio emission is detected in 
M87, where = GM/c^ is the black hole’s gravitational ra¬ 
dius (Doeleman et al. 2012). Thus, particle acceleration must 
also be occurring quite close to the black hole, where these 
processes either do not occur or are strongly suppressed. 

Within the context of M87, there are numerous spectral and 
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Name 

Eq. 

Us 

5 

iig 

18 


28 

iic 

6 

£ 

15 

ic 

10 

A 

17 

TmaxJC 

8 

Tniax,C 

11 

Tmax 

3 

777 

27 

^00 

29 

r^F 

2 

B 

44 


Table 1 

Key Symbol Definitions 

Definition 

Density of seed photons 

Number density of charges in the gap 

Number density of charges at the end of the cascade 

Cooling length of electrons due to inverse Compton scattering off seed photons 
Mean free path of 7 -ray photons due to parr production on seed photons 
Cooling length of electrons due to curvature radiation 
Gap thickness 

Maximum electron Lorentz factor due to the acceleration over the distance fic 

Maximum electron Lorentz factor due to the acceleration over the distance Iq 

Maximum electron Lorentz factor due to the acceleration over the distance A 

Lorentz factor of electrons below which they do not encounter enough seed photons to pair produce 

Lorentz factor of electrons below which they no longer efficiently Compton cool 

Lield line angular rotation frequency 

Magnetic field strength at the stagnation surface, r ~ lOr^ 


morphological constraints that should inform any attempt to 
model the lepton loading and acceleration in jets from low 
luminosity AGN (LLAGN; potentially comprising the major¬ 
ity of the population). Most importantly, at a wavelength of 
7mm and shorter the luminosity from the extended jet consid¬ 
erably exceeds that from the core, and thus importantly the lu¬ 
minosity is not core dominated. This morphology is inconsis¬ 
tent with the core-dominated morphology produced by mod¬ 
els in which the bulk of the accelerated leptons are produced 
at the black hole, e.g., models in which nonthermal energy in 
electron-positron pairs is Ug± oc (Moscibrodzka et al. 2009; 
Shcherbakov et al. 2012; Dexter et al. 2012; Moscibrodzka 
& Lalcke 2013) or with the annihilation of high-energy pho¬ 
tons produced in low-luminosity radiatively-inefficient accre¬ 
tion hows (Moscibrodzka et al. 2011). It is, however, consistent 
with models in which the energetic lepton population is gener¬ 
ated at a height of roughly lOrg (Broderick & Loeb 2009). 

This height is suspiciously similar to the location of the 
stagnation surface: the location at which the electromagnet- 
ically driven outfiow is arrested by the gravity of the black 
hole. At present, the stagnation surface is not consistently 
treated in simulations, which often assume a sufficient number 
of charged particles to maintain the electromagnetic currents 
required by the electrodynamics. However, the evacuation of 
this region (above by outfiow, below by accretion) should re¬ 
sult in the creation of large electric fields, potentially capable of 
both loading the jet with leptons and accelerating these leptons 
to high energies. 

Acceleration of particles at gaps in black hole magneto¬ 
sphere, in analogy with particle acceleration in pulsar mag¬ 
netosphere gaps, have been discussed by a number of authors 
(see, e.g. Beskin et al. 1992; Hirotani & Okamoto 1998; Levin¬ 
son & Rieger 2011), including due to field-spin misalignment 
(Neronov & Aharonian 2007). However, such studies are often 
either focused upon regions unassociated with the stagnation 
surface, require energetic backgrounds, or ignore the post-gap 
evolution of the resulting energetic particle population (Levin¬ 
son & Rieger 2011), which we find critical to reproducing the 
observed particle densities in M87. Existing efforts to numeri¬ 
cally model the post-gap process have not yet self-consistently 


modeled the gap itself (Vincent & Lebohec 2010). Thus, these 
are not immediately applicable beyond the immediate scope of 
their derivation, and none are directly relevant for the relativis¬ 
tic lepton loading of M87. 

Here we show that for M87 this is indeed expected to be 
the case. Unscreened electric fields are capable of initiating an 
inverse-Comptonpair-creating cascade. Within this process the 
number density of nonthermal particles is a function primar¬ 
ily of the soft-photon background produced by the surrounding 
accretion how and within the jet itself. Most importantly, we 
show that for M87 the inferred nonthermal lepton number den¬ 
sities match those implied by recent mm-VLBI observations 
that resolve the jet launching region. In §2 we discuss physical 
processes responsible for particle acceleration and cooling in 
the gap near the stagnation surface. In §3, we discuss the prop¬ 
erties of post-gap cascade and describe our numerical method. 
In §4 and §5 we discuss our results for M87 and Sgr A*, re¬ 
spectively. In §7 we conclude. 

2. PARTICLE CREATION AND ACCELERATION IN THE STAGNATION 

SUREACE 

The stagnation surface is continually being evacuated as ma¬ 
terial either accelerates outward in the jet or inward into the 
black hole. In the absence of some mechanism to continually 
replenish the density of charges the result is a charge-starved 
region, incapable of screening the large electric fields that will 
be present. The typical electric field strength within the stagna¬ 
tion surface is set by the angular velocity of the magnetic field 
line footprints, Up, dragged about by the spinning black hole 
or the surrounding accretion flow, 

Ef^flpRE/c, ( 1 ) 

where R is the cylindrical radius of the stagnation surface. Lor 
black hole driven jets H/r is set, in turn, by the angular velocity 
of the black hole horizon, 

HH = ac/2rH, (2) 

where eh = rg(l -E VT-t?) is the size of the horizon in Boyer- 
Lindquist coordinates. Typically, analytical studies and GR- 
MGHD simulations find rip Ri 0.25flH-0.5HH (Blandford & 
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Znajek 1977; Komissarov 2004; Tchekhovskoy et al. 2010; 
McKinney et al. 2012; Penna et al. 2013; Lasota et al. 2014). 

This unscreened electric field will necessarily accelerate 
electrons across the charge-starved gap. That is, given a gap 
thickness of A, in principle, the maximum electron Lorentz 
factor is 

7tnax~e£'A. (3) 

Eor gap sizes of order the gravitational radius, in the case of 
M87 this produces Lorentz factors of 10*^, more than sufficient 
to produce an inverse-Compton cascade. In practice, 7tnax is 
limited by radiative losses and the size of the gap is set by the 
cascade. Here we estimate the gap thickness and corresponding 
pair densities injected into the jet by the gap. 


2.1. Radiative Losses 


2.1.1. Inverse Compton 

The source of seed photons for inverse-Compton scattering is 
not completely clear. Nevertheless, we might assume that these 
are related to the observed emission from the source. Thus, if 
the observed luminosity is synonymous with the seed photon 
luminosity we will have 

Ls ^ nr^u,c/3 = LID-Fs , (4) 


where is the size of the emission region, O is the solid angle 
of the emission region (i.e., beaming factor), is the local en¬ 
ergy density of seed photons, D is the source distance, and T) is 
the observed seed-photon flux. Thus, in terms of the observed 
flux. 



3L, 


(5) 


where is the seed photon luminosity. The cooling length to 
inverse-Compton scattering is related to Us and 7niax via 


iic = 


3meC^ 

4crrMi7max’ 


( 6 ) 


where (jj is the Thomson cross section^. This defines the typi¬ 
cal distance over which particles can accelerate within the gap, 
thus 

7max,ICWreC^ Ri , (7) 


which implies a limiting Lorentz factor of 


7max,IC 


f3eRCtFB 

\ 40-7-MjC 


1/2 

= 4.2 X , (8) 


after inserting the expression for E in which we have adopted 
the standard notation, R = cm, s“\ 

5 = ^210^ G, and Us = o erg cm“^. The associated cooling 
length is = 2.2 x cm. 


2.1.2. Cun’ature 

Curvature radiation losses provide a second, seed-photon in¬ 
dependent, limit upon 7niax- The power emitted depends upon 
the curvature radius of the magnetic field, ~ R, and is given 
by (Levinson & Rieger 2011) 

p _ ^ /max /Qx 

€ 3^2 

^ This assumes that 27 ^ 3 ^^^ ~ 7 maxwhich we will verify is approxi¬ 
mately true for M87 in what follows. If this is not true, £ic is larger, producing 
a weaker limit upon 7max- 


Outflowing 



Figure 1. Schematic of a small portion of the stagnation surface and the asso¬ 
ciated gap region. Within the gap particles are rapidly accelerated to asymp¬ 
totic Lorentz factors, limited by inverse Compton cooling on the ambient soft- 
photon background (red photons), inducing a net charge separation. The result¬ 
ing up-scattered gamma rays (blue photons) pair produce on the soft-photon 
background, resulting in a pair catastrophe. When the density in Equation 
(18) is reached, particle acceleration ceases, and the remaining post-gap cas¬ 
cade produces an asymptotic lepton density at fixed charge density. Particles 
propagating downward are accreted by the black hole, while those propagating 
upward populate the black hole jet. Due to the high Lorentz factors typically 
involved, the entire process is essentially one-dimensional; the oblique propa¬ 
gation angles of the resulting gamma rays are exaggerated for clarity. Relevant 
locations and characteristic length scales are labeled (see the main text for de¬ 
tails). 


The corresponding cooling length is 

3R~meC^ 

= .,2 3 - 

2e7max 

implying a limit of 


7max,C “ 


f3R^CLFB\ 

V 2 ec ) 


1/4 

= 3.2x 


( 10 ) 


( 11 ) 
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In practice, the maximum particle Lorentz factor is set by the 
minimum radiative loss limit. As we will see in Section 4 the 
limit arising due to inverse-Compton losses is typically more 
severe. This is true if 


u, > 


r 3e^n fB\ 
\ Scr^Rc / 


1/2 

= 1.8 X erg cm“^ 




and thus just above the SED peak for our fiducial parameters. 
Hence, the inverse-Compton cooling of the gap-accelerated 
pairs initiates a pair-production catastrophe. 

The pair-production cross section peaks near 2eth at ~ 
O.SScTj (Gould & Schreder 1967). Thus, the mean free path of 
the up-scattered seed photons is approximately 


''77 


2 eth 

^5',th^77 


i 6 . 6 x 10 '°R,y^H 


R 

‘F,-4^2 


u j ? 7 thcm, (15) 


where lelO'® cm and D = D26l0“® cm. That is, as long 

as the seed photon density is sufficiently high, the gap dynam¬ 
ics is set by inverse-Compton losses. Even for the highly sub- 
Eddington sources (e.g., M87 and Sgr A*) this is often satis¬ 
fied. 


2.1.3. Synchrotron 

Generally the electric field in the gap will not be aligned 
with the magnetic field. As a consequence, accelerated par¬ 
ticles will gyrate, and thus emit synchrotron radiation. This 
can provide an alternative source of hard photons oriented at 
substantial angles to the magnetic field. The cooling length for 
synchrotron is is = iicuJuB, where ug = B~/Ett, and thus we 
generally have is <C iic, since the seed photon energy density is 
typically much smaller than the magnetic energy density in the 
sources of interest. However, note that synchrotron is only ca¬ 
pable of cooling the transverse motion of the particles. While 
this serves to ensure in extremely short time that the leptons 
are traveling along magnetic field lines (making the particle 
acceleration essentially a one-dimensional process), it does not 
otherwise limit the particle Lorentz factors. 

2.1.4. Collective Processes 


Einally, we consider the importance of collective effects due 
to the counter-streaming electron-positron plasmas. Due to 
their opposite charges, electrons and positrons will be accel¬ 
erated in opposite directions, producing a pair of cold, counter- 
propagating lepton beams. This is unstable to the well-known 
two-stream instability, modified slightly due to the opposite 
charges in the two beams (which we refer to as the “counter¬ 
streaming instability” and presented briefly in Appendix A). 
The instability growth rate is generally dependent upon the 
stream densities, corresponding to a cooling length of 


^csi 


V Ttehi^ ) 


(13) 


and thus we cannot fully evaluate the importance of plasma 
instabilities without an estimate of the pair densities. However, 
for the estimate presented below (Equation 18), ^csi ~ 7.1 x 
cm, and therefore this is never relevant. 


2.2. Pair Density within the Gap 

As the limiting radiative process, inverse-Compton cooling 
also provides a limit upon the gap thickness. The up-scattered 
photons have energies of roughly 

Thus, at meV energies, comparable to the photon energy 
near the peak of the spectral energy distribution (SED) of the 
sources of primary interest here (see Sections 4.1.4 and 5.1.3), 
the up-scattered photon energy is similar to that of the incident 
electron. Thus, the up-scattered photons will pair produce off 
of photons with energies above 


AnipC^ 


Cth ' 


: meV, 


(14) 


where th = Mi,o/?7 th is the energy density of seed photons with 
energies above eth. Eor a source with infrared spectral index a 
(i.e., Tv oc in M87 a ~ 1.2), the number of seed photons 
above an energy is oc e““- Thus, if there is a spectral break 
at Ch, below which the seed photon energy density may be ne¬ 
glected. 


f 1 eth < Ci 

l(eth/efc)“ eth>eh. 


(16) 


Eor sufficiently high e*, it is possible for p* ^ E implying that 
i-y-y ^ic- When p* = E this is larger than i^c by roughly a 
factor of three, implying in this case that the gap scale height 
comparable to either. A more quantitative model of the gap lep¬ 
ton population may be obtained in one-dimension by assuming 
that the pairs are accelerated instantaneously to their asymp¬ 
totic velocities. As shown in Appendix B.l, this implies a gap 
thickness of 

A ^/2Ec^ 7 ~ 1-7 X cm, 

(17) 

consistent with this. Generally, iy.y is larger than i\c and there¬ 
fore the asymptotic Lorentz factor is gained on a length scale 
i\c A thereby justifying our assumption that the pairs are 
instantaneously accelerated. 

Counter-propagating positrons or electrons constantly reini¬ 
tiate the cascade, which continues until sufficient charges are 
produced to screen the electric fields. Thus we would antic¬ 
ipate that the gap can be evacuated only over a single scale 
height. The top of the gap is then roughly defined when the 
charge density is sufficient to generate the gap electric held 
gradients, and thus the lepton density in the lab frame can be 
expressed as (see Appendix B.l); 


V-E 

47re 


E 

47reA 




cm 


(18) 


Note that this exceeds the Goldreich-Julian density by a multi¬ 
plicity factor given by 



(19) 


resulting from the smaller typical scale set by the gap thick¬ 
ness.^ Here we have assumed that only a single species is 
present at the top of the gap, an assumption that is justified 
by the rapid acceleration of electrons and positrons in opposite 
directions. However, the nonlinear response of the magneto¬ 
sphere will supply the charges of the opposite sign to screen 
the gap electric field and ensure that the microphysical current 
matches the global magnetospheric current (which is generi- 
cally of order of Goldreich-Julian current, jcj = enac). 

We note that first-principles particle-in-cell (PIC) simula¬ 
tions are needed in order to simulate the cascade in detail and 
properly compute the flow of charges and connection between 


^ This factor was not included in previous works (e.g., Levinson & Rieger 
2011; Moscibrodzka et al. 2011 ). 


7max 
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the pair cascade and global MHD currents. Eortunately, our 
simplified, steady-state consideration of the gap structure ap¬ 
pears to capture the most crucial aspects of the gap structure 
to within an order of magnitude. Eor instance, in PIC simula¬ 
tions of polar cascades in pulsar magnetosphere (TimokJiin & 
Arons 2013, hereafter, TA13), the simulated pair multiplicity 
Tj ~ 100 (for their polar gap thickness of A ^ 1.7 x 10"^ cm 
and a characteristic length scale R = 10^ cm; see, e.g., Eig- 
ure 22 in TA13) agrees with our estimate (19), which gives 
Tj = R/A Ri 60. Eurthermore, in agreement with the findings of 
TA13, our gap thickness is much larger than the skin depth. 


^ S V in5„l/4fi.l/4oV4 r1/4 -1/4 
^ -5.6x 10 Pjjj /(jg U^Q , 


with the skin depth given by 


dg — 


4TTe^n„ 


Encouraged by this agreement, we will adopt the results of 
our steady state gap model and consider astrophysical conse¬ 
quences. 

The relativistic outflow of particles at the top of the gap can 
carry a substantial total kinetic luminosity. Given the above 
estimates for the number density and typical Lorentz factor, 
the kinetic flux in the outflowing leptons is 

T^lep “ 7max.ICl^4gC fig 

U " 7 7 -1/2 , , (20) 

= 3.4 X 10^^Ri5H^_4B2?7th ^ erg cm ^ s ^. 

Integrating this across the entire stagnation surface gives a total 
kinetic luminosity of 

^cascade ~ f 47rR'dR'Fi^j, 

Jo (21) 

= 1.1 X erg s“'. 


Note that typically oc M ' and R cx M, where M is the mass 
of the black hole. Hence, at fixed magnetic field strength the 
pair luminosity scales as M^, implying that L^kp /Ledd oc M. Our 
fiducial numbers correspond roughly to M87, and thus a black 
hole mass of order 10*° Mq, for which the above luminosity 
is about 3 x 10“^LEdd- This suppresses the cascade luminosity 
for less massive black holes, as we will see on an example of 
Sgr A* in §5. 


2.3. Gap Stability 

The lepton structure in the gap is generally quite unstable, 
both due to the large-scale GRMHD processes that lead to its 
formation and the particle acceleration that fill it. In the case 
of its generic structure, the typical variability timescale is of 
order the light crossing time of the radial position of the gap, 
typically Ri lOrg. The radial position of the gap varies by about 
< 10% on this time scale, as can be seen in Eigure 2. 

The particle acceleration processes will induce local variabil¬ 
ity on much shorter timescales, comparable to the light cross¬ 
ing time across the gap (see Appendix B.2). In our fiducial 

model, this is roughly s. EorM87, this 

is roughly six orders of magnitude smaller than GM/c^l How¬ 
ever, in practice, the gap has a transverse extent that greatly ex¬ 
ceeds A, and therefore usually many independently fluctuating 
regions are visible at once. Assuming that each independent 



Figure 2. Meridional (x—z) snapshots of the stagnation surface, at which the 
radial velocity vanishes {u^ = 0), shown at different times with different colors 
(see legend) for a simulation of magnetically-an'ested disk around a spinning 
black hole with dimensionless spin a = 0.9 (Tchekhovskoy et al., in prepara¬ 
tion). For clarity, we only show the highly magnetized region, pmag/pc^ > 10. 
The black hole event horizon is shown as the black filled circle. Panel (a) shows 
a series of snapshots of the stagnation surface taken at intervals comparable to 
the light crossing time, Ar 5rg/c. The stagnation surface is clearly variable 
on such a short timescale, with its radial position changing by 5-10%. On aver¬ 
age, the stagnation surface is at a distance r 5-lOrg. Panel (b) shows how the 
shape of the stagnation surface evolves on longer time scales, Ar lO'^rg/c. 
On such longer time scales the amplitude of variability can be as high as 50%. 
For example, the shape of the surface can be substantially distorted by infalling 
gas that pushes the jet (e.g., see cyan dash-dotted line, which corresponds to 
t = 75305rg/c). 


region has a local radius of A, if the visible extent of the stag¬ 
nation surface is Ry, the net degree of variability on the particle 
variation timescale is roughly 


5ng A 
(«g) Rv 


( 22 ) 


Eor low frequencies, Ry ~ R, and this reduction is enormous. 
However, for the very-high energy gamma rays produced di¬ 
rectly by inverse Compton scattering, the intrinsic relativistic 
beaming will substantially reduce Ry. Eor emission associ¬ 
ated with electrons with Lorentz factor 7 , a patch of the surface 
that will be visible has a typical size of R/ 7 , implying that for 
7 > 10 ° the observed particle density variations, and therefore 
the flux variations, will be of order unity on the particle accel¬ 
eration timescale. In our fiducial model, where the soft seed 
photon energy density peaks around 1 meV, this corresponds 
to a gamma-ray energy above a few GeV. 

As the energy of a gamma-ray increases, the number den¬ 
sity of soft photons that lie above the pair creation threshold 
increases, and so does the optical depth to pair produc¬ 
tion. The existence of the pair catastrophe implies that at suf¬ 
ficiently high e.y the optical depth to pair production is high. 
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i.e., ^ 1. Since the energy threshold for annihilation of a 

gamma ray is cx , the number density of soft photons above 
the pair creation threshold (and therefore the optical depth to 
annihilation) scales as cx e“. The normalization of the total 
gamma-ray optical depth can be obtained in terms of set 
by = 1 for gamma rays with energy 7 max,ic«JcC“: 


r {c )- ^ ^ 4OT^c‘*/e7) 

^77 ^ 4n7gC^/7maxJc) 

Hence, for gamma-rays with energies 


r ( e/OTgC^ \ “ 

^77 \ Tmax.IC / 

(23) 


7 ^ ^7,min — Tmax.IC^rZ^C 



(24) 


the seed photon distribution will be optically thin to the result¬ 
ing inverse Compton gamma rays, allowing the latter to escape. 
Eor values typical of M87, this corresponds to energies of or¬ 
der 67 ~ 630 GeV, above which significant attenuation can be 
expected. 


3. PARTICLE CREATION AND ACCELERATION BEYOND THE 
STAGNATION SUREACE 

The leptons streaming out of the gap opened near the stag¬ 
nation surface will still have near-maximal Lorentz factors, 
Tmax.ic- As a result, they will continue to pair produce beyond 
the gap each time they traverse a distance of roughly iic- How¬ 
ever, unlike the results of the previous section, since the leptons 
are outside of the gap, they are no longer rapidly re-accelerated, 
and thus the pair catastrophe is no longer being driven. The 
result is a redistribution of the energy of gap-injected leptons 
among a set of daughter leptons produced via the resulting final 
cascade. 


3.1. Qualitative Asymptotic Estimates 

We have already noted that for gamma rays with energies be¬ 
low a lower limit determined by the seed photon SED the ambi¬ 
ent soft photon population will be transparent, halting the pair 
production cascade. Explicitly, for a power law seed photon 
distribution above a low-energy cutoff e* (as described above 
Equation 16) the fraction of gamma rays with energies above 
£ 7 ,min following a single scattering is 


/>£.,,mi„ ( 7 ) 


(e7,n,in/27^ei) “ 7 < ,min/2£;p (25) 

1 Otherwise, 


where we have assumed that e* is sufficiently small. This is 
necessarily a strong function of the initial lepton Lorentz factor, 
dropping rapidly for Lorentz factors below \J and 
thus e 7 _mm implies a limit upon 7 . 

However, since the energies of the upscattered gamma rays 
are typically small in comparison to the lepton energy near this 
limit, i.e., 27 ^ei -C ymc^, in practice the lepton will undergo 
many scatterings before cooling appreciably. That is, there will 
be roughly 


Ascat 


jmc^ mc^ 
2y^eb 2yeb ’ 


(26) 


opportunities to produce a gamma ray with energy above £ 7 ,min. 
This increases with decreasing 7 and e*, ameliorating the 
decrease in />, due to the latter. Thus the number of 
pair producing gamma rays is expected to be approximately 
Ascat/>e.y min ■ Setting this number to unity provides the desired 


limit upon the lepton Lorentz factor, below which electrons 
cease to typically have sufficient energy to continue the pair 
cascade: 


•yl-OL A-a.^ot 


l/( 2 a-l) 


777 


(27) 


Eor our fiducial model, 777 ~ 7.7 x 10"^. This limit on the re¬ 
distribution of energy implies an associated limit on increase 
in the particle density. An even redistribution, maximizing the 
number of leptons produced, results in 

IIoo 7max,IC /oo\ 

- , 

rig 7 -^-y 


and is roughly 6000 daughter leptons for each lepton produced 
by the gap for our fiducial model. 

The absence of the generation of energetic gamma rays does 
not imply the absence of additional inverse-Compton cooling. 
Hence, a second, typically more stringent limit arises from 
inverse-Compton cooling directly. Eor the cooling length to 
become comparable to the scale height of the seed photon dis¬ 
tribution requires 


3meC^ - 

-> r 

4aTMs7oo 


< 

^ Tmax,IC 

r 




(29) 


We note that in our fiducial model, the cooling of the leptons 
from to is unable to produce gamma rays energetic 
enough to pair produce. Therefore, the resulting radiation can 
escape and can be a substantial source of emission. 

Since yoonoo ^ 7 max,icng most of the energy imparted to the 
leptons by the accelerating electric field across the gap has been 
radiated via the Comptonization of the soft seed photons. Typ¬ 
ical energies range from the optical to the GeV, with an SED 
that depends upon that of the underlying soft photon popula¬ 
tion. However, this emission will be highly beamed along the 
magnetic field lines that govern the particle motions. 


3.2. Numerical Asymptotic Estimates 

In practice the estimates obtained in the previous section 
depend critically upon assumptions regarding the final elec¬ 
tron spectrum, that is itself dependent upon the soft pho¬ 
ton SED. Eor this reason we also computed the asymptotic 
number and energy distribution of electrons produced due to 
inverse-Compton scattering off of a given seed photon SED 
from a nearly mono-energetic injection at high Lorentz factor. 
Throughout this section, for compactness we will measure en¬ 
ergies in units of the electron rest mass energy. 

We make two simplifying assumptions, both of which are al¬ 
most certainly well justified. The first is that we may ignore 
any subsequent pair annihilation due the comparatively high 
lepton energies and low lepton densities. The second is that the 
number density of soft seed photons vastly exceeds the num¬ 
ber densities of up-scattered gamma rays, a consequence of the 
short mean free path to pair production and the low lepton den¬ 
sities. These remove the non-linear terms in the coupled Boltz¬ 
mann equations for the leptons and the gamma-rays. Thus, it 
is sufficient to consider the evolution of a single lepton through 
multiple Compton generations. That is, we compute the energy 
probability distribution of a single lepton as a function of the 
number of scatterings instead of time. The asymptotic form 
of this probability distribution then describes the asymptotic 
spectrum of a full population of leptons. 

Explicit inputs are the seed photon spectrum and electron 
injection energy, which we leave as a parameter, eg. Eor the 
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former we assume a power law model between some minimum 
(Cm) and maximum (e^) energies, consistent with the observed 
radio and infrared emission from the sources of interest, i.e.. 


dN, 

Me) = ^ oc 
de 



ejn ^ e Cm 
otherwise. 


(30) 


the threshold e.y,min, above which the seed photon distribution 
is optic ally-thin to pair production. 

In principle, we should then set cm by the pair production 
optical depth; from Equation (24) this gives a typical 

Cm ~ 4m2cV(630 GeV) (35) 


The values of Cm and cm are determined by the particular source 
under consideration, though typically when a large enough dy¬ 
namic range is chosen the resulting evolution becomes insen¬ 
sitive to their specific choice.^ The linearity in the gamma-ray 
distribution is synonymous with assuming that remains fixed 
throughout the post-gap pair cascade. 

With the above we may now compute the update in the lep¬ 
ton energy distribution with each scattering. We do this in two 
steps, first computing the implied up-scattered gamma-ray dis¬ 
tribution associated with a mono-energetic electron distribu¬ 
tion, and then from that inferring the changes to the electron 
distribution arising from inverse-Compton losses and pair pro¬ 
duction. To compute the former, we begin with the relation¬ 
ship between the original seed photon energy Cs, up-scattered 
gamma-ray energy e-^, and initial electron energy Cgf. 

Cs{c^:Ce,i)=\'^eej{ce^i-C^) ^ ’ (31) 

I 0 otherwise. 


which for our fiducial parameters gives cm~ L6 eV. In prac¬ 
tice, we find the simulation results are insensitive the particular 
value we choose for cm- In all of our simulations, unless spec¬ 
ified otherwise, we choose cm = 0.8 eV (as we discuss below, 
results for cm = L6 eV are essentially identical; see Eigure 6). 
In a similar fashion to Equation (34), we can write down the 
distribution function of photons that are below the pair produc¬ 
tion threshold, c-y < e-y,rnin; these photons escape the system 
without undergoing scattering and thus form the high-energy 
emission spectrum, 

pOO 

fy'!yj+Ml)=‘^ / dCe,igic^,Ce,i)ff'^iCe,i)Q{c^.mm-e.y). (36) 
Jo 

Combining the above, we obtain the energy distribution of 
the leptons and photons in the (j+ l)th generation via 

f^Me) = fl^jM^e)+f!;iM^e) , (37) 

= (38) 


where we assumed Cs Cej and approximated energy of the 
electron post-collision as Cej = Cej-Cy. Equation (31) may be 
obtained from the standard Compton scattering formula in the 
high-energy electron and oblique seed photon limits. As a re¬ 
sult, for a monoenergetic electron distribution, the up-scattered 
gamma-ray distribution is 


gie^jCe^i)— fslCsiCyjCej)] — 


/j[c.s(C7 1 £«,<)] 

2(Ce,i-cM ■ 


(32) 


Beginning with some initial injection distribution, we can com¬ 
pute the asymptotic distribution, in which generations corre¬ 
spond in a loose sense to time or height from the injection 
point. 

Eor this, we integrate Equations (37)-(38) numerically. We 
start the integration at the top of the gap, and inject a single 
electron at an energy 

Eq — 'Tmax.ICn^eC , (39) 


As mentioned above, since <C Cej,Cj, we approximate the 
electron energy following scattering by Cej ~ Ce.j-c^. There¬ 
fore, if fj{ce) is the electron distribution after the y'th scatter, the 
electron distribution is rearranged by inverse-Compton losses 
to 

pOO 

flCj+l(^e)= / deejg(eej-ee,eej)ff'^(eej). (33) 
Jo 

In practice the integral need only be computed up to the initial 
maximum energy of the electron distribution. 

Similarly, since Cj^2^ Cs, i.e., much greater than the rest 
mass of the produced pair, which is itself much greater than that 
of the seed photons, we may approximate the energy of each 
lepton in the resulting pair by Cg = c^jl. Hence, the distribution 
of new pairs is 

poo 

M (ee) = 4 / dc,,ig(2ce, CeMf^’i^e,) , 

Jo 

(34) 

where 0(.r) is the Heaviside function and ensures that the in¬ 
tegral is only over regions in which the gamma-ray can pair 
produce on the soft background. Note that we assume that all 
gamma-rays will pair produce as long as their energy exceeds 

^ The reason is that when /tQ the up-scattered photons will be 

unable to pair produce (and the energy loss becomes comparatively insignifi¬ 
cant) and when cm 3 > 4mjC^/eo the number density of photons with energies 
near em will be small in comparison to those responsible for the bulk of the 
particle creation (and thus may be neglected). 


which is the energy gained by an electron between its conse¬ 
quent encounters of seed photons in the gap and is therefore the 
characteristic energy with which the electrons emerge from the 
gap. We discretize Equations (37)-(38) on a logarithmic grid 
in energy, which extends from Emin = 10 “^OT(,c^ to Eniax = 2£o- 
Eor numerical reasons, we smooth out the initial energy distri¬ 
bution by 1.5% so it is properly represented on our numerical 
grid, by choosing the following initial injection distribution. 


fo\e,)^ 


dCe 


1 

\/27rcr| 


exp 


jCe-Eo)^ ' 

2al J’ 


(40) 


with (T£ = 0.015£o- 

In order to ensure robust evolution of pair cascade numeri¬ 
cally, we found it to be important to ensure lepton number con¬ 
servation under Inverse Compton cooling (Equation (33)). To 
do this, we compute the discretization of the integral in Equa¬ 
tion (33) on two energy grids that contain the same number of 
grid cells but are shifted relative to each other by half a cell. 
If the changes in the number of leptons due to these two dis¬ 
cretizations are different by more than 50% relative to each 
other, then we use their weighted sum that preserves the total 
number of leptons exactly. Otherwise, we use the discretization 
that results in the smallest change in lepton number. We imple¬ 
mented the numerical code in a Cython-based Python library. 
In order to speed up the integrations, we parallelized the code 
using OpenMP. We describe the numerical results of cascade 
evolution in Section 4.3. 
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4. APPLICATION TO M87 

Thus far we have remained agnostic regarding the particular 
sources for which particle acceleration at the stagnation surface 
is relevant. Here we apply this model to the generation of the 
nonthermal particles in the radio jet of M87. We begin with 
a short overview of the relevant source parameters and their 
corresponding gap parameters. 

4.1. M87 Source Parameters 

Observations of stellar dynamics within M87’s sphere of in¬ 
fluence produce an inferred mass of 6.6 x IQ^Mq, assuming 
a distance of 17.9Mpc (Gebhardt et al. 2011).'* The energy 
output of M87 peaks at 1 mm (300 GHz or 1 meV), at which 
point the flux at Earth is roughly 1 Jy. The typical isotropic- 
equivalent total luminosity is L r; 10 "*^ erg s“* Ri 10“®LEdd, 
where LEdd is the Eddington luminosity. This makes M87 quite 
underluminous and dominated by emission at wavelengths 
about the peak frequency, as we discuss below. The typical 
scale near the black hole is the gravitational radius of the black 
hole, Pg ~ 10'^ cm. Here we discuss the empirical implications 
of various additional observations for the physical properties 
near the stagnation surface. 

4.1.1. Jet Velocity and Orientation 

Apparent motions near the black hole of roughly 2c have 
been observed (e.g.. Walker et al. 2008). These imply a min¬ 
imum Lorentz factor of E r; 2.3 by distances of 0.3 mas, cor¬ 
responding to linear distances of 8 x 10 *® cm r; 80rg, and im¬ 
plying an inclination of 25°. On scales of 1", more than three 
orders of magnitude farther out, apparent motions of 6c have 
been observed, implying E > 6 at an inclination of r; 10°. Here 
we will employ the radio implied Lorentz factor and associated 
orientation, E = 2.3 and 0 = 25°, where specific values are re¬ 
quired. 

4.1.2. Magnetic Field Strength Estimates 

The observational implications for the magnetic field in M87 
depend upon the structure and dynamics of the material in the 
vicinity of the black hole. Here we present two estimates, as¬ 
suming in the first that the observed mm-wavelength emission 
is due to a radiatively inefficient accretion flow (RIAE, Narayan 
& Yi 1994, 1995), and in the second that the jet is Poynting 
dominated and responsible for the power observed in the ra¬ 
dio lobes. Both reach similar conclusions regarding the held 
strength. 

As suggested by their name, RIAEs are characterized by their 
extraordinarily low radiative efficiency, i.e., L = with 

Prf R! 10“^-10“"*. Since most of the gravitational binding energy 
is released, and therefore most of the luminosity is produced, 
near the black hole, this implies an accretion flow density of 


47rrg2/3c Aitridrg^ 


thus we arrive at the first magnetic held strength estimate of 

1 /2 

B r; r; ( r; 3 x G. 

V Bdfg C / 

(42) 

It is not clear that the mm-wavelength emission is associ¬ 
ated with an accretion flow. At 7 mm the jet dominates the 
source morphology (Walker et al. 2008). Thus, we present 
a second estimate based upon the assumption that the jet is 
Poynting dominated near the black hole and relating the me¬ 
chanical power at small and large radii. The Poynting lumi¬ 
nosity associated with black-hole driven jets is (Tchekhovskoy 
et al. 2010): 

= (43) 

where k k. 0.045 for parabolic jet geometry, 4> is the poloidal 
magnetic flux threading the black hole, and Hh is the angular 
velocity of the horizon of the black hole, which encodes the 
dependence upon the black hole spin.® 

A variety of estimates of the total jet power in M87 have been 
obtained across a wide range of scales. Radio lobe measure¬ 
ments probe emission extending over 30 kpc, and have resulted 
in estimates of Ljet Ri 10"**ergs“' (de Gasperin et al. 2012). This 
is similar to estimates arising from the modeling of the complex 
of knot features on kpc-scales (Owen et al. 2000), and on 60 pc 
scales from interpreting the HST-1 complex as a recollimation 
shock (Stawarz et al. 2006). Thus, we adopt this estimate as a 
conservative estimate of the total jet power. Eor a characteris¬ 
tic black hole spin, a = 0.9, equating these then gives a poloidal 
magnetic field strength estimate at the event horizon: 

Bp = y r; 2.5 X 10^ G, (44) 

where A = 27r(a^r^-l-3r^)/3 is the area of the one of the two 

hemispheres of the event horizon, where th = r'g(l -I- V 1-a^) is 
the radius of the black hole event horizon. 

Despite the significantly different physics being invoked, the 
two estimates are quite similar. While this is not accidental 
(thick disks are believed to be a necessary component of jet 
formation), it does critically depend upon the low radiative ef¬ 
ficiency generally implicated in vastly sub-Eddington sources 
(such as M87). 

At the distances of the stagnation surface, which lies inside 
the light cylinder, the magnetic held is poloidally-dominated. 
If we now assume that the jet is roughly parabolic, i.e., Rj/rg r; 
(r/rg)'/^, within the gap region the above estimates imply 

Br.^r.35G. (45) 

for a typical stagnation surface distance of r ~ lOrg. 

4.1.3. Synchrotron Cooling Times 

In the jet frame, the observed synchrotron cooling timescale 
is 


where we have assumed a disk height of r and an inflow veloc¬ 
ity near the horizon of jSc. The jet magnetic held will accumu¬ 
late until it balances the ram pressure of the accretion flow, and 

^ While quantitative differences appear if we were to adopt the lower mass 
of 3.5 X 10® Mq implied by gas dynamical modeling (Walsh et al. 2013), the 
qualitative consequences remain unchanged. 


, 3mc^ 6nmc^ j Vg 3 / 2'Kemc \ 

= 4aTu'gCj' " V 

^ Note that the low-spin approximation given by Equation (43) remains ac¬ 
curate for spins a < 0.95, beyond which it over-estimates the true power by 
up to 30%. For a higher order expansion that remains accurate for all spins, 
see Tchekhovskoy et al. (2010). See also Tchekhovskoy et al. (2012) for the 
comparison of vaiious approximation to jet power. 
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where v' and v'g = eB' jlmneC are the observation and cy¬ 
clotron frequencies in the jet frame, and 7 ' is the Lorentz fac¬ 
tor of nonthermal particles, again measured in the jet frame. 
To relate this to emission in the lab frame, note that v' = 
^r(l-/3cos0) and since the magnetic field is likely to be 
dominantly toroidal, B' k. B/T. Hence, as measured in the lab 
frame. 


t =t r! 

'sync - 'sync^ ' 


3 

Gj 


liremcT^ 1 


B^v{\- PcosQ) 


(47) 


1.8 X 10 "Bi 


5 „- 3 / 2 ,- 1 / 2 ^ 


^ 11.5 


where B = SOBi.sG and v = SOOj^h.sGHz and the remainder of 
the quantities are given by our fiducial M87 values. 

Evident in Equation (47) is the strong dependence of fsync 
on B and E, and therefore position. As the jet propagates B 
decreases and E increases in a fashion that depends on the spe¬ 
cific jet structure, resulting in a corresponding rapid increase 
in fsync. Generally, magnetic flux conservation implies that the 
poloidal and toroidal components of the magnetic field within 
the jet scale oc R~^ and oc R~^, respectively (it is for this rea¬ 
son that despite having similar strengths near the black hole, 
jet magnetic fields quickly become strongly toroidally domi¬ 
nated). Eor a parabolic magnetic field, R oc which we as¬ 
sumed when obtaining Equation (47). Similarly, E is expected 
to rise oc R. However, in practice the Lorentz factor of the jet in 
M87 appears to asymptote to r; 5 at large radii, limiting its im¬ 
pact. Regardless, near the stagnation surface these imply that 
fsync will initially grow oc eventually slowing to ex 

The implication for the image structure is obtained by com¬ 
paring fsync to the outflow timescale. 


which grows linearly with jet height. When fsync > 4pp the rel¬ 
ativistic leptons cannot efficiently cool, and therefore the non¬ 
thermal particle density is effectively conserved. At moderate 
distances from jet base (z > ISfg) the above estimates suggest 
that this is the case. At much larger distances where the Lorentz 
factor saturates (z > lO^fg) the synchrotron cooling becomes 
efficient. 

Near the stagnation surface fsync is comparable to foutflow 
However, this is complicated by the large initial bulk Lorentz 
factors resulting from the acceleration across the gap. As a 
result, a substantial fraction of leptons still resides in the or¬ 
dered non-thermal relativistic component that emerged from 
the polar cascade propagating along magnetic field lines and 
therefore geometrically immune to synchrotron cooling. The 
study of how these nonthermal electrons isotropize and cool is 
beyond the scope of this paper. Thus, for simplicity we will 
ignore synchrotron cooling in the estimation of the energetic 
lepton population of the jet. 


Lutflow ■ 


/3c 


3.3 X 10= 


4.1.4. Seed Photon Density Estimates 

The source of the soft seed photons in M87 is not immedi¬ 
ately clear due to potential contamination from the jet itself. 
Nonetheless, it is now certain that the millimeter-wavelength 
emission is associated with compact features with scales com¬ 
parable to that of the horizon (Doeleman et al. 2012). This pro¬ 
vides circumstantial evidence that the observed millimeter, in¬ 
frared, and optical emission arise very near the black hole itself. 
Thus, even if the emission is from within the jet itself, at the 


jet heights of relevance the bulk Lorentz factors are expected to 
be small (of order unity), and hence not strongly beamed. As a 
consequence, we estimate of the seed photon density from the 
observed emission directly. 



V (Hz) 


Figure 3. Spectral energy distribution for M87 from the radio to the opti¬ 
cal, constructed from the literature (Spencer & Junor 1986; Baath et al. 1992; 
Junor & Biretta 1995; Sparks et al. 1996; Lonsdale et al. 1998; Biretta et al. 
2002; Ly et al. 2004; Kovalev et al. 2005; Lister & Homan 2005; Petrov et al. 
2007; Perlman et al. 2007; Lee et al. 2008: Baes et al. 2010; Doeleman et al. 
2012, additional details will be published elsewhere). Green filled triangles 
show a subset of the infrared-optical flux limits collected in Broderick et al. 
(2015). Filled circles and open triangles con'espond to flux measurements in 
which the source is resolved and unresolved, respectively; we take the filled 
circles as directly indicative of the near-horizon seed photon distribution and 
the open triangles as upper limits. Errorbars indicate variability, not intrinsic 
measurement uncertainty. For comparison, a simple broken power-law with 
high-frequency spectral index a = 1.2, as employed in the text, is shown for 
reference. 


The SED of M87, shown in Eigure 3, is reasonably well fit 
by a broken power law, i.e., oc where a evolves rapidly 
as the source becomes optically thin (see Eig. 3): 


fO < 300 GHz 
[1.2 z/> 300GHz. 


(49) 


The flux density at the break is roughly 1 Jy. If all of the emis¬ 
sion above 300 GHz arises in a radiatively inefficient accretion 
disk, then the soft photon energy density is roughly 

Us Ri 0.01 erg cm“^, (50) 

where we have now fixed the location of the stagnation surface 
to lOrg. This is well above the densities necessary for inverse 
Compton losses to dominate the limits upon the acceleration of 
the leptons within the gap. 


4.1.5. Implied Nonthermal Particle Properties 

If the sub-mm luminosity arises due to optically thin emis¬ 
sion of the jet launching region, as suggested by the recent 
mm-VLBI observations (Doeleman et al. 2012), it implies a 
rather soft nonthermal particle density, with dn/d^ k, 7 “^-^. 
The source size is roughly 40 pas, implying that the associated 
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millimeter flux is related to the particle density by 


Stt^c 

= 4.9x 10“2 


V 

f—V 

\40/ras ) 


'■ 2aj^ 


2a. 
min p 


a +1 


lOr, 


a 11 
2~^~6 




a 1 
—I— 
2 6 


Qrn 


(51) 


where vb = eB/27rmeC is the cyclotron frequency, r(x) is the 
standard gamma-function, i/ = 300 GHz, and we have assumed 
a typical lower cutoff on the accelerated particle distribution 
of 7 min ~ 10^. Comparing this to an observed flux of roughly 
1 Jy and assuming the magnetic field estimate in Equation (45) 
gives an approximate nonthermal particle density estimate of 


n Ri 15 cm' 


(52) 


consistent with efforts to quantitatively model the emission re¬ 
gion of M87 on horizon scales (Broderick et al., in prepara¬ 
tion). 


4.2. Particle Acceleration at the Stagnation Surface 

We can now estimate the properties of the gap, and the as¬ 
sociated accelerated particles, from M87’s parameters given 
above. This gives our fiducial model for M87. The asymp¬ 
totic lepton Lorentz factor is determined by inverse Compton 
losses and is roughly 


7max — 7max,IC ~ 1 .6 X 10 , (53) 

and this sets our fiducial value of electron energy, Eq. Whereas 
its numerical value is known only to a factor of few, we choose 
not to round this value in order to get a self-consistent answer. 
This corresponds to a seed photon pair-production threshold of 

eth ~ 1.2 meV A* ~ 1.1 mm. (54) 

Therefore, photons at and just below the spectral break in 
M87’s SED will contribute to pair production within the gap. 

The length scale over which the asymptotic lepton Lorentz 
factors are obtained is (see Equation 6 ) 

^icRi5.5x 10'°cm = 5.6x IQ-Vg, (55) 

and the typical mean free path of the up-scattered photons to 
pair production is (see Equation 15) 

1.7x lO'^cmRi 1.7x 10“Vg. (56) 

Hence, as anticipated, both the acceleration and subsequent 
pair production within the gap occurs on scales much smaller 
than the typical gap thickness that might be expected on global 
considerations (i.e., r^). The resulting gap thickness is then 
roughly (see Equation 17) 

A R! 4.3 X 10" cm, (57) 

implying a corresponding density at the top of the gap of (see 
Equation 18) 

Ug R 2.2 X 10“^ cm“^. (58) 

We note that since A by about an order of magnitude exceeds 
£ic, over which the pairs attain their terminal Lorentz factor, 
our earlier assumption that the pairs instantaneously attain their 
terminal Lorentz factor is justified. 


4.3. Post-Gap Cascade 

Given the threshold seed photon energy of 1.2 meV, the min¬ 
imum gamma-ray energy for which the seed photons bath is 
optically-thick to pair production is e-y min = 640 GeV, hence 
the asymptotic Lorentz factor for M87 is (see Equation 27) 

7 .y.y = 2.6x 10®, (59) 

implying enhancement in the number density due to the post¬ 
gap cascade of roughly noajn^ = 670 (see Equation 28). 

To verify this analytical estimate, we carried out numeri¬ 
cal integration of the pair cascade equations (37)-(38), as de¬ 
scribed in Section 3.2. We choose the following parameters 
to describe the seed photon spectrum in our fiducial model of 
M87: a= 1.2, e,„ =1.2 meV, tu = 0.8 eV. We start with an ini¬ 
tial distribution that corresponds to an electron emerging from 
the top of the gap at an energy Eq. Lor numerical convenience, 
we represent its energy distribution with a narrow Gaussian 
distribution shown in Eigure 4(e)-(h) with thin solid dark red 
line. Different rows show distributions for different values of 
Eo/nieC^: 10®, 10^, 1.6 x 10® (fiducial value), and 10'®. We 
then evolve the system of equations (37)-(38) numerically, as 
described in Section 3.2, and show subsequent distributions in 
Eigure 4 with colored lines (see the legend). 

We adopt as our fiducial model pair cascade for M87 the 
case when Eq = 1.6 x lO^m^c^, shown in Eigure 4(g). The 
electron cools (see Equation 33) due to the upscattering of a 
seed photon into a 7 -ray photon, which in turn pair produces 
on an additional seed photon (see Equation 34). These two 
contributions to lepton energy distribution are clearly visible in 
Eigure 4(h): after the first scattering, i.e., in generation 1, the 
IC-cooled electron distribution is a power-law extending from 
7iep ~ 0.5£ 7 ,min = 6 X 10® to 7 iep = 2 X 10® and the new pairs 
form a Gaussian distribution centered around 7 iep r 7 x 10®. 
Note that the lower-energy cutoff in the power-law emerges be¬ 
cause the seed photon bath is transparent to gamma-rays with 
energy less than e-y,min- 

Numerically-obtained values of density enhancement vs. 
generation number are shown in Eigure 5(a), for different val¬ 
ues of Eq. It is clear that the density enhancement saturates 
around generation Agen ^ few x 10. The numerical result for 
density enhancement in our fiducial model of M87, n^jng = 
640, is shown with long-dashed magenta line and is in excel¬ 
lent agreement with the analytic expectation, discussed above 
and given by Equation (28), «oo/ng ~ 670. This yields a jet 
lepton density near the stagnation surface of 

tioo ~ 15 cm“® . (60) 

In the absence of nonlinear plasma phenomena, these energic 
pairs will inverse-Compton cool over roughly a single jet scale 
height, i.e., within a distance r, at which point they will have 
undergone large-angle deflections. These leptons, now mis¬ 
aligned with the background magnetic field, will subsequently 
synchrotron cool and isotropize. The rate of this cooling de¬ 
pends on the microscopic Lorentz factors and becomes compa¬ 
rable to the outflow timescale only when 7 r 10 ““® on scales 
of the gap height. Thus the resulting population is in excellent 
agreement with that required to produce the observed sub-mm 
emission. 

As is also clear from Eigure 5(a), higher initial energy leads 
to a higher density enhancement. Similarly, Eigure 5(b) shows 
that the total energy and the energy in leptons are higher for 
higher Eq- The flatness of the solid lines in Eigure 5(b) demon- 
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Figure 4. A sequence of the resulting energy distributions for escaping photons (left column) and leptons (right column) for different values of initial energy of the 
electron, shown in different rows, from top to bottom: E{)/meC^ = 10^, 10^, 1.6 x 10^, 10^^. Different generations of lepton and photon energy spectra are shown with 
lines of different colors, thicknesses and line types (see the legend). We start at generation zero with a single electron at energy Eo. The electron up-scatters a seed 
photon into gamma rays and cools into a power-law energy distribution. The gamma ray pair-produces off of another seed photon (see main text for details). Gamma 
rays, which are not energetic enough to pair-produce, escape, and their spectrum is shown in the left column. Note that at generations of Pi 100 and greater, the shape 
of the spectrum of leptons and photons is essentially independent of the initial energy of the electron, E^. it is only the normalization of the spectra that depends on 
£” 0 , not their shape. 


strates that the total energy of the cascade (carried by leptons 
and photons) is conserved at high accuracy by our numerical 
scheme. Note that as seen in Eigure 4, there exist asymptotic 
properties of the pair cascade that are independent of Eq. Eor 
example, it is explicitly apparent in Eigure 5(c) that the av¬ 
erage lepton Lorentz factor, (7iep), is independent of Eq for 
Agen ^100 and is controlled completely by 

We investigated the dependence of density enhancement fac¬ 
tor on key parameters of our model, and the results are shown 
in Eigure 6. The blue dots connected by solid lines show the 
numerical results; they are in excellent agreement with the 
analytical scaling given by Equation (28). The density en¬ 
hancement is most sensitive to the initial electron energy, Eq-. 
Aiep,oo oc Eq [see Eigure 6(a)], and the gamma-ray energy below 
which the seed photon bath is optically-thin to pair production, 

£ 7 ,min: Aiep.oo OC ~ (^^6 Eigure 6(c)). However, 

if we now recall that both Eq = 7max,ic«reC^ and e-y,niin oc 7max,ic 
depend on 7max.ic (see Equation (24)), we find that these de¬ 
pendencies nearly cancel out, and that density enhancement is 
most sensitive to the density of seed photons or source lu¬ 


minosity: 


n„ 


( \(a-l)/(2a-l) i/(2o!-1) ^ ( \0.14 

OC (7maxJCefcj nj^ ' Ri (7max.ICeij « 


0.71 


( 61 ) 

Importantly, the numerical results are insensitive to the upper 
energy cutoff of the seed photon spectrum, as shown in Eig¬ 
ure 6(d). 

We ensured that our simulation results are numerically con¬ 
verged. Eor this, we varied the number of energy bins for more 
than 3 orders of magnitude and computed the error relative to 
the highest-resolution simulation carried out using 2x10® en¬ 
ergy bins. As seen in Eigure 7, the simulation results converge 
at first order with increasing energy resolution. Eor the simula¬ 
tions described here, we employed a fiducial resolution of 10^ 
energy bins, for which the relative error is less than 0.3 per cent 
in all quantities. 

The asymptotic Lorentz factor after the pair cascade satu¬ 
rates, set by when the inverse-Compton cooling length is com¬ 
parable to the jet scale height, is roughly 10^. This implies 
direct synchrotron emission up to « 10-100 eV, well into the 
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Figure 5. Evolution of various quantities in pair cascade with generation for 
different values of initial lepton energy, £ 0 ^ which is measured in units of mgC^. 
Panel (a) shows the dependence on the generation number of density enhance¬ 
ment, A^iep = f^co/fig- Panel (b) shows the dependence of the total energy in 
leptons and photons with colored solid lines and the lepton total energy, ^lep, 
with other lines, for different values of initial lepton energy, Eq. Panel (c) 
shows the evolution with the generation number of the average lepton Lorentz 
factor, ( 7 iep>. 


ultraviolet. 


4.4. Direct Inverse Compton Signal 

The direct inverse Compton signal from the stagnation sur¬ 
face itself should be visible up to cvhegr ~ 640 GeV, above 
which it will be absorbed by the pair cascade. As is clear from 
Eigure 4, the direct IC signal from the cascade has a flat energy 
spectrum in E^dN/dE, with every decade in energy carrying 
approximately the same amount of energy. This is consistent 
with the analytical expectation given in Appendix C. Thus, it is 
possible that the observed very-high energy gamma-ray emis¬ 
sion is associated with the stagnation surface (see, e.g., Levin¬ 
son & Rieger 2011). Emission in the band pass of FermiPLPIE, 
around e.^ ~ 10 GeV, is mostly due to Compton up-scattering of 
seed photons near the peak of the SED, Cj « 1.2 meV, by elec¬ 
trons with Lorentz factors 7 fa (e.y/2e,5)'/2 = 2 x 10®, which is 
comparable to 7 .^-^. The associated power per lepton is 

Rio GeV ~ /!>. (62) 



Figure 6. Dependence of density enhancement factor, noo/f^g, on the parame¬ 
ters of the cascade. Connected blue dots show simulation results and the green 
dotted lines show the analytical expression given by Equation (28). Panel (a) 
shows the dependence of density enhancement factor, «oo /ng on the initial en¬ 
ergy Eq of the lepton, panel (b) on the lower energy cutoff in the seed photon 
distribution €m, panel (c) on the minimum energy of gamma-rays able to pair 
produce e-y ^in (see Equation 24), and panel (d) on the upper energy cutoff in 
the seed photon distribution. The agreement between the model and the simu¬ 
lation is better than 30%. 



Figure 7. Dependence of fractional error on resolution in our fiducial model 
(evaluated at 40th lepton generation) for various quantities: lepton number, 
A^iep, emitted photon number, A^ph, lepton energy, E’lep, emitted photon energy, 
Eph. The error is computed relative to a high-resolution model with 2 x 10^ 
spectral energy bins, and exhibits linear convergence. At our fiducial resolution 
of 10‘^ energy bins, the relative en'or is less than 0.3 per cent in all quantities. 


The highly relativistic nature of the lepton distribution pro¬ 
duces a large relativistic aberration, beaming the up-scattered 
gamma rays within a cone with an opening angle 1 around 
the original electron momentum vector, and thus concentrating 
the emission in a solid angle of However, it also 

limits the region of the stagnation surface that can be viewed 
to an area Finally, inverse-Compton cooling also sets a 

length scale over which we expect to find particles that can pro¬ 
duce sufficiently energetic up-scattered gamma rays, ^ic.io Cev 
typically much smaller than rg, and here roughly 0.05rg (see 
Equation 6 ). Thus, the typical anticipated 10 GeV flux pro¬ 
duced near the stagnation surface is independent of the beam¬ 
ing, assuming that the viewing angle is within the jet opening 
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angle: 


e 110 GeV 


10 Gev 


Vl^R] 


ttRI 


‘oo ^2 ^jcao GeV = 


(63) 


= 7.1 X 10 erg cm ^ s '. 


In principle, all leptons with Lorentz factors above 7 will con¬ 
tribute to the energy flux. In practice, this increases the flux 
near 10 GeV by a factor of less than two (see Appendix C), to 




12(a-l) , ttR] 

10 Gev- 2a-1 

= 1.2 X 10“^ erg cm“^ s“'. 


(64) 


This is considerably larger than the (7.5 ±5.9) x 
10“'^ erg cm“^ s“' reported for M87 by Fermi in the first 
FermHV^P^ catalog of > 10 GeV sources (lEHL, Ackermann 
et al. 2013). It is, however, consistent with the distance- 
adjusted fluxes of other gamma-ray bright blazars in the Fermi 
sample (see, e.g., Ackermann et al. 2011). 

The disparity is almost certainly due to the high degree of 
beaming anticipated. At 10 GeV, the high energy emission 
is beamed within 0.1” of the original electron direction. Be¬ 
yond this angle, the energy of the up-scattered photons drops 
dramatically; at 10 ° the typical up-scattered photon energy is 
only roughly 70 times larger than that of the original seed pho¬ 
ton, well below the energies of interest. Thus, the gamma- 
ray emission from the stagnation surface would be expected to 
drop precipitously outside of the jet half-opening angle. This 
is a strong function of height, scaling cx near the black 
hole, though at our fiducial height of the stagnation surface 
this is 0j = R/z — 18°. Because the jet is collimating the tan¬ 
gent to the field lines, which correspond to the range of angles 
over which the relativistic leptons are directed and subsequent 
gamma rays emitted, makes a considerably smaller angle with 
the jet axis, — 0.59j = 9° for a parabolic jet. Both of these 
are well within the inclination implied by radio observations 
of 25°, suggesting that it in this case it would be rare to find 
substantial gamma-ray emission in M87. 


4.5. Total power of the cascade 

Using Equation (21), we obtain the total power dissipated by 
the cascade, 

^cascade - lO"^^ ergs“', (65) 


corresponding to roughly 10% of the total jet power, Ljet — 
10'*'* ergs“*. Thus, most of the jet power flows out in the form 
of Poynting flux and only 10% of it gets converted into parti¬ 
cle energy at the stagnation surface and the cascade. Most of 
this power is radiated in the form of gamma-rays and is not 
observable due to beaming away from our line of sight, as we 
discussed above. A small fraction of this power. 


Too^oo J _ Too J 

^cascade “ ^cascade 

Tmax,IC^^^ 

i 4 X 10“^Lcascade = 5 X 10"*° ergs“*, 


( 66 ) 


is carried asymptotically by the cascade and will be radiated via 
synchrotron emission. This is roughly a factor of 20 smaller 
than the isotropic-equivalent luminosity, consistent with the 
beaming correction inferred from detailed source modeling 
Broderick Sc Loeb (2009).® This is not surprising since the 


® Note that since the angular size of the source is fixed in the estimates of 
the seed photon densities, the beaming conection does not affect Hj. 


asymptotic values of density and Lorentz factor produced by 
the cascade, Uoo and 700 , are consistent with those necessary to 
explain M87’s radio emission, as discussed in § 6 . 

We note that in addition to total energy budget our model has 
the potential to address an important morphological feature of 
M87 jet that is pronounced in radio images at a frequency of 
43 GHz (Walker et al. 2008): edge-brightening of the jet. The 
electric field in the gap, given by Equation (1), increases away 
from the rotational axis. This leads to density enhancement 
peaked toward the edge of the jet (see Equations 8 and 61) and 
therefore has the potential to lead to stronger jet emission near 
the jet edges. More detailed study of this issue is warranted. 

5. APPLICATION TO SGR A* 

A second obvious application is to the putative jet in Sgr A*, 
the black hole at the center of the Milky Way. Unlike M87, 
there is no obvious radio jet in Sgr A*, making the task of con¬ 
straining its properties somewhat less well defined. However, 
there is some recent evidence for a collimated outflow contain¬ 
ing a non-relativistic population of leptons (Li et al. 2013). In 
this case the putative jet lights up on parsec scales as a re¬ 
sult of intersecting, and shock heating upon, a stream of gas 
in Sgr A*’s vicinity. Again, we will first review some relevant 
parameters for Sgr A* and then assess the implications of the 
particle acceleration model presented in Section 2. 

This is complicated by the current uncertainty regarding the 
dynamics and morphology of the region responsible for the ob¬ 
served sub-mm emission. This uncertainty manifests itself pri¬ 
marily in the assumed magnetic field strength. In sections 5.1- 
5.3 we will assess in detail the implications for a putative jet as¬ 
suming that the observed emission arises primarily in a RIAE. 
In section 5.5 we consider the possibility of a self-consistent 
solution in which the jet itself is responsible for the sub-mm 
emission. 


5. 1. Sgr A * Source Parameters 

The mass of and distance to Sgr A* are the best known of 
any black hole candidate, obtained from observations of the 
orbits of individual stars, yielding M = 4.3 ± 0.4 x 10 ®Mq 
and 8.3 ±0.4 kpc, respectively (Ghez et al. 2008; Gillessen 
et al. 2009). As with M87, the SED peaks near 1 mm, 
where the typical observed flux is roughly 3 Jy. However, 
the isotropic-equivalent luminosity is much lower, roughly 
L K, 10^® erg s“* Ri lO'^Ledd- Hence, even in Eddington units, 
Sgr A* is considerably more underluminous than M87, for 
which L R! lO'*^ ergs”* r 10“®LEdd- Here we collect various 
parameter estimates relevant for the putative jet launching re¬ 
gion near Sgr A*. 

We will find that Sgr A* and M87 have roughly similar pho¬ 
ton energy densities, the frequency of the spectral peak, and 
magnetic field strengths. This leads to comparable (to within a 
factor of few) physical (in cm) gap thicknesses in both cases. 
However, the mass of Sgr A* black hole is about 1500 smaller 
than in M87. This means that in SgrA*, the gap thickness is 
comparable to large-scale size of the system, and this leads to 
the suppression of the cascade. 

5.1.1. Jet Velocity and Orientation 

In the absence of an obvious jet, there are no direct measure¬ 
ments of a putative jet velocity. In contrast, if the jet feature 
reported in Li et al. (2013) is real, it is necessarily oriented at 
a large angle relative to the line of sight. Here we will assume 
a similar, moderate Lorentz factor for the jet (E r 2.3) and an 
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inclination angle of 90°, where necessary, consistent with the 
assumed parameters in Li et al. (2013). 


5.1.2. Magnetic Field Strength Estimates & Synchrotron Cooling 

Times 


The detection of polarized emission at mm-wavelengths, and 
thus the measurement of the Earaday rotation measure of « 6 x 
10^ radm“^, produces an estimate for the accretion rate of M = 
2 X 1 O“®M 0 yr“*-2 x 10“’M© yr“\ depending on the details of 
the accretion model assumed (Agol 2000; Marrone et al. 2007). 
Eollowing equations (41)-(42) gives a corresponding estimate 
of the possible confined jet magnetic field of up to 


B 



10’“^G. 


(67) 


In both cases, the accretion power, ranging from 
10 ^^- 10 ^°ergs“', exceeds the potential jet power by two 
orders of magnitude, and imply the need for low radiative effi¬ 
ciencies, i.e., rid ~ 0.01-10“"^. These field strengths are similar 
to those inferred from spectral fitting (e.g.. Yuan et al. 2003; 
Broderick Sc Loeb 2006). Eor concreteness we will assume 
the lower field strength, corresponding to the comparatively 
higher radiative efficiency, i.e., we will assume 100 G on the 
horizon and thus roughly 10 G near the stagnation surface. As 
a consequence, the synchrotron cooling timescale is similar to 
that of M87: 

fsync « 4.0 X s. ( 68 ) 


Since the gravitational radius of Sgr A* is much smaller, so is 
the apparent outflow timescale, comparable to 


Gutflow 


2.4 X lO’ 



s, 


(69) 


and thus for Sgr A* synchrotron cooling is completely irrele¬ 
vant to the evolution nonthermal particle population on horizon 
scales. 


5.1.3. Seed Photon Density Estimates 

Unlike M87, Sgr A* exhibits a clear sub-mm bump, of¬ 
ten modeled as a thermal disk component (Yuan et al. 2003), 
though also potentially produced by the jet component asso¬ 
ciated with the jet launching region (Ealcke & Markoff 2000; 
Yuan et al. 2002; Moscibrodzka & Ealcke 2013). Nevertheless, 
above 1 mm the SED is well modeled by a power law with 
a ~ 1.25 (see Eigure 8 ); 


r 0 ^ < 300 GHz 

\l.25 j/> 300GHz. 


(70) 


Eurthermore, short-timescale variability from the mm through 
the infrared imply that this emission arises near the black hole, 
and thus may be identified with the soft seed photons relevant 
for particle production and acceleration at the putative jet stag¬ 
nation surface. 

The implied seed photon energy density above 1 mm is 

Us ^ 1.5 X 10“’ erg cm“’, (71) 

which is comparable to that implied in M87 (see Equation 50). 
We will And, however, that the relevant seed photon population 
lies above 10*’ Hz (wavelengths shorter than 0.26 mm), for 
which the energy density is reduced by roughly a factor of 5. 



Figure 8 . Spectral energy distribution for Sagittarius A* from the radio to the 
optical, collected from Yuan et al. (2004, and references therein) and Marrone 
(2006). Filled circles and open triangles con'espond to flux measurements in 
which the source is resolved and unresolved, respectively; we take the filled 
circles as directly indicative of the neai'-horizon seed photon distribution and 
the open triangles as upper limits. Errorbars indicate variability, not intrinsic 
measurement uncertainty. For comparison, a simple broken power-law with 
high-frequency spectral index a = 1.25, as employed in the text, is shown for 
reference. 

As we will see below, this leads to about a factor of 3 increase 
in the physical thickness of the gap (in cm) in comparison to 
M87. This turns out to lead a qualitative change in the cas¬ 
cade operation: since Sgr A* black hole mass is ^ 1500 times 
smaller than in M87, Sgr A* gap thickness comes out to be 
comparable to the size of black hole, and this suppresses the 
cascade. 


5.1.4. Implied Nonthermal Particle Properties 

The roughly power-law, declining SED of Sgr A* in the sub- 
mm again justifies assuming the emission is well characterized 
by optically thin emission from a soft population of nonther¬ 
mal electrons. In this case, dn/d^ oc 7 “’'^, quite similar to 
M87. With a source size of 37 ^as, equating the expression in 
Equation (51) to the observed millimeter flux of 3 Jy yields a 
nonthermal particle density estimate of 

na!l0®cm“’. (72) 

This is much higher than for M87 as a direct result of the much 
smaller mass, and therefore more compact emission region. 

Given the previous magnetic held estimate, the typical 
Lorentz factors are expected to be of order lO’. While the 
infrared emission requires considerably larger 7 ’s, it is also 
highly variable and thus potentially associated with additional 
dissipative events. 

5.2. Particle Acceleration at the Stagnation Surface 

As might be expected given the similarities in the soft photon 
densities, the inverse-Compton limited Lorentz factor is about 
a factor of 2 smaller than that found in M87 (see Equation 3; 
see Equation 53 for M87 value): 

7max~7.6x 10'. (73) 

The length scale over which this is attained is (see Equation 6 ) 
fic~8.4x 10 '° cmRi0.13rg. 


(74) 
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This value in cm is larger by about a factor of 2 and in units 
of gravitational radius by about a factor of 2300 than the M87 
value (see Equation 55). As a consequence of the moderately 
lower 7niax the threshold energy for seed photons which are 
able to pair produce on the subsequently up-scattered inverse- 
Compton gamma is moderately higher, 

eth = 2.7 meV A* ~ 0.46 mm, (75) 

which is pushing into the far-infrared. Hence, unlike in M87, 
for Sgr A* the soft photons responsible for pair production in 
the gap are significantly above the spectral break. 

Due to the larger eth, the corresponding mean free path to 
pair production is comparatively large (see Equation 15), 

K. 6.6 X 10*2 cm Ri lOrg . (76) 

This value exceeds that for M87 by a factor of ~ 4 in cm and 
in units of by a factor of 5900 (see Equation 56). This im¬ 
plies that at most, the gap pair catastrophe is limited to a single 
generation. This suggests that the global jet dynamics, will 
play a significant role in the determining the structure of the 
gap. Ignoring this potential complication, the equilibrium gap 
thickness is of order (see Equation 17; see Equation 57 forM87 
value) 

1.1 X lO^^cmft! L7rg, (77) 

itself larger than a gravitational radius (though given the wide 
disparity between £ic and this may be highly variable). The 
associated density at the top of the gap is (see Equation 18) 

ng Ri 2.5 X 10“^ cm“^, (78) 

about an order of magnitude smaller than for M87 (Equa¬ 
tion 58). 

In summary, the gap densities for Sgr A* (Equation 78) and 
M87 (Equation 58) are within an order of magnitude of each 
other because of the similarity of seed photon densities, mag¬ 
netic field strengths, and the peak frequency of the spectrum. 
The order of magnitude smaller density in Sgr A* comes from 
a factor of few smaller values of magnetic field strength and 
above threshold photon density for Sgr A*, falling well short 
of the «oo — 10^ cm“^ inferred from radio observations (see 
Sec. 5.1.4). 

53. Post-Gap Cascade 

Since is much larger than the system scale height at the 
gap, i.e., i-f-y ^ r Ki lOrg, there is effectively no post-gap cas¬ 
cade. This is because for compact emission regions the pair- 
production optical depth falls significantly at large radii for two 
reasons. Eirst, the seed photon density drops oc r~^ for suf¬ 
ficiently large r. Second, the direction of propagation of the 
seed photons and gamma rays become increasingly aligned as 
they travel far from their origin, reducing their center-of-mass 
energy and increasing the energy threshold for pair production 
oc r-. Hence a gamma ray that fails to annihilate in the first 
mean free path is unlikely to ever do so. 

Thus, iioo ~ Ug ~ 2.5 X 10“^ cm“^, implying that in this case 
it is likely that Sgr A* emission is likely coming from the ac¬ 
cretion disk or from the disk outflow just outside of the highly 
magnetized jets (e.g., Moscibrodzka & Ealcke 2013). If a sub¬ 
stantial fraction of Sgr A* arises from jet emission (but note 
that this does not have to be the case: in fact, SgrA* might not 
have any jet at all), some other particle loading and accelera¬ 
tion mechanisms must be occurring, such as gas entrainment 


from the accretion disk and the acceleration of electrons via 
reconnection and/or shocks. 

5.4. Gamma-Ray Emission 

According to Equations (20) and (21), the pair cascade lu¬ 
minosity scales as oc M^ng7niax.ic- Since black hole mass is 
smaller by 3 orders of magnitude, density is smaller by an or¬ 
der of magnitude, and 7tnax,ic is smaller by a factor of two, 
the pair cascade luminosity is roughly ~ 7 orders of magni¬ 
tude smaller than in M87 (Equation 65), or Ri 3 x 10^^ ergs“'. 
Since 7oo ~ 10^ the overwhelming majority of this is deposited 
in gamma rays beamed along the jet with typical energies of 
7max»Lc^ ~ 400 TeV. 

Based on the overall gamma-ray luminosity it is tempting 
to relate this to the power in the 1 GeV-100 GeV gamma-ray 
10 kpc-scale jet features found by Su & Einkbeiner (2012). 
Such an identification, however, requires at least two addi¬ 
tional mechanisms. Eirst, the gamma rays must be reprocessed 
down in energy by between three and five orders of magni¬ 
tude. Second, the emission must be scattered in a distributed 
fashion along the jet to produce the observed extended feature. 
Accordingly, we leave the discussion of potential signatures 
of gamma-ray emission from Sgr A*’s putative jet for future 
work. 

5.5. Self-Consistent Jet Solutions 

We finish our discussion of Sgr A* with a consideration 
of potentially self-consistent jet solutions. While it has been 
found above that a jet confined by an observationally motivated 
RIAE models is expected to produce little intrinsic emission, 
this need not be the case if the RIAE picture is relaxed. 

Increasing the magnetic field increases the efficiency of the 
inverse-Compton pair catastrophe. Eor B > 200 G becomes 
small in comparison to the gap height, and thus multiple gen¬ 
erations of pairs can be sustained resulting in a well-defined 
inverse-Compton pair catastrophe. Simultaneously, the number 
of pairs required to reproduce the observed sub-mm emission 
decreases. EorB r; 10"^G the two are approximately equal, sug¬ 
gesting the existence of a self-consistent jet model similar to 
that found for M87. As with weaker fields synchrotron cooling 
may continue to be ignored (this is true for B < 1.5 x 10"^G). 

However, note that 10"^ G would be much larger than even 
optimistic estimates for the magnetic field strength in Sgr A*. 
The implied jet power would be roughly Lem ^ 10"^*ergs“*, 
five orders of magnitude larger than Sgr A*’s bolometric lumi¬ 
nosity. While such large luminosities have been implicated in 
Sgr A*’s past (e.g., the Eermi bubbles Su et al. 2010), it is dif¬ 
ficult to imagine sufficiently low radiative efficiencies to hide 
such jet power in the present epoch. 

Perhaps more problematic is the need to confine such a large 
magnetic flux. This is ostensibly done via an accretion flow, 
which implies a mass accretion rate of 

B^r ^ 

MR^>10-^MQyr-', (79) 

comparable to the Bondi rate for Sgr A*. As discussed in Sec¬ 
tion 5.1.2, such high accretion rates are already excluded by po- 
larimetric observations. Moreover, the corresponding radiative 
efficiency of the accretion flow would need to be exceedingly 
small: < 10“®, much lower than the already small efficien¬ 

cies implicated in RIAE models. 

Eor these reasons we consider a self-consistent jet model 
driven by a gap-powered inverse-Compton pair catastrophe to 
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be highly disfavored as an explanation of Sgr A*’s current 
sub-mm emission. Nevertheless, this may not be the case for 
Sgr A*’s recent past, where the luminosity of Sgr A* has been 
inferred to be considerably higher. 

6 . DISCUSSION 
6.1. Source of seed photons 

Our model relies on the presence of seed photons with en¬ 
ergies e ,5 Ri 1 meV, or wavelengths A ~ 1 mm, to sustain the 
cascade. The spectra of Sgr A* and M87 above the peak of 
the SED give us an estimate of the seed photon density at the 
source. However, so far we have been agnostic about the ori¬ 
gin of these photons; whereas in the case of Sgr A* they could 
plausibly only come from the accretion disk, in M87 they could 
come from the disk, the jet, or both (see, e.g., Broderick & 
Loeb 2009). 

In our model of M87, the synchrotron luminosity of the jet, 
which is lit up by the leptons accelerated in the cascade, natu¬ 
rally explains the SED of M87 above its peak at Cj > 1 meV. 
This is because the pair cascade catastrophe naturally leads to 
pair density of Hoo ~ 15 cm“^ and Lorentz factor ~ 10"^ 
(see §4.3 and §4.5): both of these values are comparable to 
the those required to account for M87’s radio emission above 
the peak of the SED (Broderick &. Loeb 2009). This means 
that in principle no external source of seed photons is required: 
the cascade is self-sustaining, similar to polar cascades in pul¬ 
sar magnetospheres (see, e.g., Timokhin & Arons 2013). That 
is, given a stray seed photon (e.g., from an accretion disk), the 
cascade fires up and produces energetic electrons that then cool 
and produce sub-mm photons that can serve as seed photons for 
the cascade. 

6.2. Comparison to other work 

An alternative process of pair production in jets is annihila¬ 
tion of gamma-rays from the disk (Levinson & Rieger 2011; 
Moscibrodzka et al. 2011). In this scenario, a sub-mm pho¬ 
ton, emitted by the disk, undergoes a sequence inverse Comp¬ 
ton scattering events with disk’s own hot electrons. These 
interactions upscatter the photon to gamma-ray energies, and 
eventually above pair production threshold. Collisions of these 
gamma-rays in the funnel mass-load the jet with pairs. 

In order for this process to lead to an interestingly large num¬ 
ber density of pairs (larger than the Goldreich-Julian density, 
ugj), the following two conditions must be met: (i) the mm- 
emission of the disk must be compact and marginally optically 
thick to Thompson scattering; (ii) disk electrons must be very 
hot, so that on average no more than Ri 1.5 Compton scatter¬ 
ings are needed to reach pair-production threshold. To satisfy 
these requirements, disk radiative efficiency is assumed to be 
high, r]d Ri 0.3 (Moscibrodzka et al. 2011). 

In contrast, our model allows for standard interpretation of a 
radiatively inefficient accretion flow in M87, <g; 1, which we 

believe is a more natural assumption. The seed photons around 
the SED peak are produced in the jet, outside of the disk proper, 
and thus they have a lower probability of being up-scattered to 
gamma-ray energies by disk electrons. Both of these factors 
combined render ineffective the process of jet pair loading due 
to 7-7 collisions in the context of our model. 

Levinson & Rieger (2011) considered pair cascade at the 
stagnation surface of black hole magnetosphere and focused 
on the gamma-ray emission from the resulting leptons. In our 
work we take a broader approach and construct quantitative ID 
numerical simulations of the post-cascade evolution. We con¬ 


sider both the direct gamma-ray signal as well the synchrotron 
emission of the leptons produced in the cascade. In fact, we 
argue that the seed photons are produced in the cascade itself 
and not the accretion flow. We use the observed spectrum of 
M87 to determine the seed photon spectrum and eliminate the 
sensitivity of the model to the uncertain details of the accretion 
flow. 

7. CONCLUSIONS 

Particle creation at the stagnation surface of Poynting- 
dominated jets presents a natural mechanism for Ailing black 
hole jets with nonthermal particles near the horizon, as re¬ 
quired by recent mm-VLBI observations of M87. This 
model provides an excellent quantitative description of a po¬ 
lar electron-positron cascade that can All black hole jets with 
relativistically-hot plasma with sufficient abundance to explain 
the observed jet emission in M87. 

Within the stagnation surface of astrophysical jets, large 
electric fields are anticipated, arising due to the charge star¬ 
vation within the jet. These are capable of accelerating stray 
leptons to extremely high Lorentz factors, limited by inverse- 
Compton cooling on the ambient soft-photon background, due 
to local synchrotron emission within the jet or through their 
production in an accretion disk. The up-scattered gamma rays 
are sufficiently energetic to initiate a pair catastrophe that ulti¬ 
mately self-consistently determines the structure and content of 
the jet near the stagnation surface, i.e., within the gap, beyond 
which sufficient charge densities exist to screen the accelerat¬ 
ing electric fields. Typical gap thicknesses are much smaller 
than the gravitational radius, and thus largely independent of 
the global jet structure. 

Subsequent to exiting the gap, the leptons typically have very 
large Lorentz factors, and thus participate in a post-gap inverse- 
Compton cascade. While this does not increase the charge den¬ 
sity, and therefore does not affect the jet structure, it does sig¬ 
nificantly increase the lepton density and decrease their spe¬ 
cific energy. Typical lepton density enhancements can reach 
10 ^, and thus the post-gap cascade can have a dramatic im¬ 
pact upon the particle content of astrophysical jets. The result¬ 
ing pairs are then natural candidates for the relativistic lepton 
populations that produce the synchrotron emission in the cores 
of radio AGN. In particular, particle creation at the stagnation 
surface provides an excellent quantitative explanation of the 
horizon-scale emission observed in M87. In contrast, it vastly 
under produces the lepton content of Sgr A*’s emission region, 
consistent with the lack of an obvious jet in that source. 

Gamma-ray luminosity produced by the stagnation surface 
pair cascade considerably exceeds that detected by Eermi/LAT 
for M87, though it is similar to those found for gamma-ray 
bright blazars. However, this naturally results from the strong 
beaming along the jet of the up-scattered emission, implying 
a viewing inclination larger than Ri 9°. Thus, within the con¬ 
text of the stagnation surface pair catastrophe model we present 
here, the lack of bright high-energy gamma-ray emission from 
M87 favors the larger inclination implied by high-frequency 
radio observations. Nevertheless, future work should com¬ 
pute the angular dependence of the emergent inverse-Compton 
spectra and the potentially strong constraints on jet orientation 
that may result. 

The authors thank Dimitrios Giannios, Ramesh Narayan 
and Alexander Philippov for helpful discussions. The authors 
would also like to thank the anonymous referee for many help- 


HORIZON-SCALE LEPTON ACCELERATION IN M87 


17 


ful suggestions that have resulted in a much improved presenta¬ 
tion. A.E.B. receives financial support from Perimeter Institute 
for Theoretical Physics and the Natural Sciences and Engineer¬ 
ing Research Council of Canada through a Discovery Grant. 
A.T was supported by a Princeton Center for Theoretical Sci¬ 
ence Eellowship and by NASA through the Einstein Eellowship 
Program, grant PP3-140131. The simulations presented in this 
article used computational resources supported by XSEDE al¬ 
location TG-AST100040 on NICS Kraken and Nautilus and 
TACC Lonestar, Longhorn, Ranch, and Stampede. A.T. thanks 
Perimeter Institute for hospitality and financial support of three 
visits to Perimeter Institute during which most of the work on 
this project was carried out. 


APPENDIX 

A. RELATIVISTIC COUNTER-STREAMING PAIR TWO-STREAM 
INSTABILITY 


One might be concerned with the potential for heating of the 
counter-streaming electron and positron beams, ultimately ra¬ 
diated away as synchrotron emission. Here we estimate the 
growth rate of this variant of the classic two-stream instabil¬ 
ity. Note that in this case we really do have a one-dimensional 
problem; particles are strongly confined to field lines and thus 
can only participate in collective motions in the direction of 
motion. 

The electrons and positron motions are governed by the 
Boltzmann equation, which upon some simplification gives the 
following relation for the charge perturbations in terms of the 
unperturbed lepton distribution functions (see, for example, the 
appendix of Broderick et al. 2012): 
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where in the second expression we have assumed that all the 
vectors are co-linear due to the afore mentioned dimensionality 
of the problem. Inserting = (n /2)(5^(p ± po) and performing 
the trivial integrals results in 
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Prom Gauss law (V • E = 47rp) we then find 
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This is minimized for the second root (the - sign), when 
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Thus, the growth rate of the relativistic ally counter-streaming 
pair instability is. 


Tcsi = 3(w): 
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Within the gap the density given in Equation (18) sets a char¬ 
acteristic upper limit, 

V-E 

-- - "Ss 


cm P (18) 
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Note that here the characteristic length scale is A, i.e., much 
smaller than the distance to the stagnation surface, lOr^. The 
corresponding plasma frequency is then 

ujp^ 1.0 X Hz , (A8) 


and thus setting 7 = 7 max, 

Ecsi ~ 4.2 X 10-®Mx,oWHz, (A9) 

which is typically exceedingly slow in comparison to the light¬ 
crossing time of the gap. The corresponding cooling length is 
^csi ~ c/Ecsi- 


B. ONE-DIMENSIONAL, APPROXIMATE GAP STRUCTURE 

The gap structure is set by the competition between pair pro¬ 
duction and the screening of the externally applied electric field 
by the resulting charge separation. In principle, this requires a 
solution to the fully non-linear coupled Maxwell-Boltzmann 
equations, including inverse Compton scattering. Particle-in¬ 
cell (PIC) simulations could be one such approach. However, 
in practice, a rough estimate of the gap structure may be ob¬ 
tained via an approximate form of the Maxwell-Boltzmann 
equations, in which the non-linearity is manifested by bound¬ 
ary conditions. 

Here we explicitly construct the solutions to the simple gap 
particle population model listed in Section 2.2. We assume 
that the electric field within the gap Eg is sufficient to instan¬ 
taneously drive particles to their inverse Compton limited ve¬ 
locities, and thus the particle distributions are essentially mo- 
noenergetic. The resulting particle distribution will evolve to 
screen the externally applied Eg, preventing further accelera¬ 
tion and saturating the growth of the gap lepton density. Note 
that in this way the linear model encapsulates the nonlinearity 
inherent in the full Maxwell-Boltzmann equations. 

Eor concreteness, we assume that the electric field is pos¬ 
itive in the upward direction, corresponding to the direction 
of the guiding magnetic field line, in which case the electrons 
and positrons accelerate downward and upward, respectively. 
Then, the number densities of the pairs («(.+, rig-) and upward 
and downward propagating gamma-rays (n.y+, n^-) evolve due 
to pair production and inverse Compton scattering via: 


i;j‘^-(2k^+\)i;Y-iY(\-]Y) = o 


=A- uY = k^+-±\ [k''-+-] +k\\-e-) 
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where = a~^.,ins, ctic = crj-nj, and iis is the number den¬ 
sity of seed photons. Pairs are generated by both upward and 
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downward propagating photons, though as assumed travel only 
in their prescribed directions. Thus, the upward (downward) 
propagating gamma-rays are produced by inverse Compton 
scattering by only the positrons (electrons). 

While Equation (Bl) is manifestly one-dimensional, the 
stagnation surface itself extends transversely over a macro¬ 
scopic scales, i.e., typically much longer than lie or The 
strong beaming of the up-scattered gamma rays by the highly 
relativistic, accelerated leptons, implies that regions that are 
separated by a transverse distance of more A/ 7 , where A is 
the gap height to be determined below, execute independent 
pair catastrophes. As a result, the stagnation surface is com¬ 
prised of a large number (rs A^) of essentially indepen¬ 
dently evolving one-dimensional, pair cascades, each described 
by Equation (Bl). 

The pair cascade is inherently unstable, as will be shown 
in detail in Section B.2 though is clear from the nature of the 
quenching: pairs form and accelerate until the resulting charge 
segregation screens the externally applied electric field. As a 
consequence, the microscopic gap structure, i.e., that within 
a single causally connected tube, is highly variable. How¬ 
ever, the macroscopic gap structure, comprised by very many 
causally connected tubes, is well described by an ensemble 
average over microscopic tubes at different stages of the dis¬ 
charge cycle, i.e., the time-averaged structure of a single mi¬ 
croscopic region. Due to the linearity of Equation (Bl), this 
corresponds immediately to the stationary solution, which we 
address first. 


B.l. Stationary Structure 

Assuming stationarity, all the time derivatives vanish, and 
Equation (Bl) is a set of coupled ordinary differential equa¬ 
tions for the gap structure: 


(nj+ + nj-'j 
-d^iig- = (nj+ + n-y-^ 

rTu-y+ = o.icne+ — (y.^ryn^+ 
O-U-y— - 


(B2) 


from which, with the required symmetry, we obtain the solu¬ 
tions 

g{z)=Aco^\\(yj2aica^^z\ . (B9) 


We may now construct / by integration: 

f = A, 
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The contribution to the electric field from the pairs within the 
gap may be obtained directly by integrating Gauss’ law. 


Ey = 4TTe(ne+-ne-) = 47re [f(z)-f(-z)] 
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When A = 0, i.e., the charge density vanishes, the total electric 
field is simply that due to the gap. Eg. The particle cascade 
saturates when the gap is closed, i.e., Ey r; 0 at the center, and 
Ey = Eg at the gap boundary, where it is presumably screened 
by MHD processes, taken to lie at a height Zg. The first condi¬ 
tion gives 

AireA 

C = -. (B12) 

Otic 

The second then sets A in terms of the gap size: 

, Q:icE„ 


47 re(cosh( 5 - 1) ’ 


(B13) 


where S = y^laicxr^Zg = 2zg/A, where A = is the 

characteristic gap length scale. 

Because all particle production is assumed to occur within 
the gap, only outgoing particles exist above and below the gap 
boundaries. The condition that / > 0 everywhere within the 
gap then requires that the density of each species vanishes at 
a corresponding boundary, set by their respective directions of 
motion, i.e., the positron and electron densities vanishes at the 
bottom and top of the gap, respectively. This implies. 


B = A 



sinh i5, 


(B14) 


Choosing z = 0 to lie in the gap center, symmetry requires that from which we infer the particle densities at the top of the gap. 


rieAz) = ne-{-z) = f(z) 


and similarly 


n^*{z) = n^-(-z) = g(z). 


Thus, Equations (B2) reduce to 

f = [.g(D+.g(-D] 
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It is clear from the above that the natural gap scale is A. 
However, the density at the top of the gap is a weak function of 
b, with sinh^/(cosh(5-1) dropping to 2 at (5 ~ 1 and falling be¬ 
low 1.3 for 5 >2. Thus, independent of the precise gap thick¬ 
ness, the density at the top of the gap is 


and 


g'= oticf - oi'^-ig, (B 6 ) 

where for compactness we have denoted dz with primes. Tak¬ 
ing the second derivative of the second equation and inserting 
the first, we obtain the following for g: 

g" + [g(z)+g(-z)] = 0. (B7) 

Noting that g'{-z) = -[g(-z)]' and g"{-z) = [g(-z)]", this im¬ 
plies 

[giz) +g{-z)f -2aica^.y [g(z) + g(-D] = 0 

y (Ijo) 

oi-i.^[g{z)-g{-z)] =0, 




47reA ’ 


(B16) 


which is directly comparable to Equation (18). Kinetic plasma 
simulations are desirable to verify this linearized model but are 
beyond the scope of this work. 

Equation (B16) gives a value of charge density that is a factor 
E/A = lOr^/A ^ 10^ times greater than the Goldreich-Julian 
density, pcj ~ HE/27rc. Consequently, the gap structure im¬ 
plies a large pair multiplicity that is, however, not unusual 
for magnetospheric cascades, e.g., in the context of pulsars 
(TA13). In fact, the multiplicity 77 = E/A implied by Equa¬ 
tion (B16) is in good agreement with the simulation results of 
TA13 (see the discussion after Equation 18). 
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B.2. Gap Dynamics 

The gap model described above is already linearized, sim¬ 
plifying the discussion of its dynamics considerably. Here we 
look at small, harmonic perturbations about the steady state 
configuration. We have some freedom to measure time and dis¬ 
tance in convenient units, which to simplify the following ex¬ 
pressions, we do in = 1 joLr^^c. Then, defining C! = uj/a-yyC, 
K = kja^^c, and ( = aic/ctj^, the equation for the perturba¬ 
tions becomes 


—/((55 — K)5ne+ = Shy -I- Sn^- 
+ K)Sn^+ = SuY + Sn^- 
—i(U5 — K)SnY = (^6ne+ — Sn^ 
-i((S + K)Sn^- = (^Sug- - Suj-. 


(B17) 


This may be trivially rearranged to define the modes as a stan¬ 
dard eigenmode problem: 



0 i 



^ Sne7 


^ Sne+\ 

0 

—K i 

i 


dng- 

= 03 

Sug- 

K 

0 K-i 

0 


Sfl'y+ 

Sfl'y+ 

VO 

iC 0 

-K-iJ 




\Sn^-J 


(B18) 


The associated four mode-specific dispersion relations are 



1± 


Y^1+4(-4k2±4\/C^-(1-I-4C)«:2 


(B19) 


where every combination of the two ± may be chosen. Two 
of the four modes are unstable (— and — 1 -), both of which cor¬ 
respond to the pair catastrophe (related to the symmetric and 
anti-symmetric components). The maximum growth rate is 


The post-gap lepton population has a height-dependent, 
mono-energetic distribution: 


— =n^S[y-y{z)], 

d'y 
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where 7 ( 2 ) is set by inverse Compton cooling, and determined 
by 


dy A-ajUs 2 
dz 3mc^ ^ ’ 


(C3) 


subject to the initial condition that at the top of the gap, 7 = yoo- 
While this may be trivially integrated, only the above expres¬ 
sion is required below. 

Given a mono-energetic electron distribution, the energies 
of the up-scattered gamma ray and the initial seed photon are 
related by 


e dcs 1 

^ 27 ( 7 - 6 ) de 2 ( 7 - 6 )^ ’ 

where we have expressed these in a slightly different form than 
usual, obtaining the seed photon energy in terms of the gamma- 
ray energy. 

The resulting gamma-ray flux is then 


f , , dndns de^ n 


(C5) 


where the factor of dcs/de simply converts from per unit seed 
photon energy to per unit gamma-ray energy. The beaming cor¬ 
rections (enhanced emission along the direction of the electron 
momentum and reduced viewing area) cancel as described in 
the main text. Affecting the integral over 7 is trivial due to the 
mono-energetic nature of the electron distribution, hence. 


r__|^^o = a^-,cS(G3__) = 2^--c. (B20) 

Thus, the growth time. Ell, is roughly the light crossing time 
of the open gap. The physical origin of the instability is easy 
to understand: additional leptons induce additional inverse 
Compton scattered photons which, in turn, produces additional 
electrons, causing the process to exponentiate. 

C. DIRECT INVERSE COMPTON ELUX DENSITY ESTIMATES 

Here we construct a simple model with which to estimate the 
direct inverse Compton signal from the gap-accelerated lep¬ 
tons. The critical input assumption is that the leptons evolve 
solely due to inverse Compton cooling after leaving the gap. 
Simplifying assumptions, that certainly hold in the case of M87 
and are likely to hold genetically, include a broken power-law 
distribution for the seed photons and a uniform seed photon 
density over the cooling length scales of interest. As in sec¬ 
tion 3.2, for compactness we will measure energies in units of 
electron rest mass. 

Eor concreteness, the seed photon distribution is take to be 
of the form given in Equation (30) 

^ ^ (a-l)(l-/3) J(e^/ei)“(^+'^ e, < e;, 

dcs el a-13 £,>£*, 

i.e., a broken power-law with a spectral indexes of a and /3 
above and below the break, respectively. The above is normal¬ 
ized such that the seed photon energy density is Mj. 
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is the electron Lorentz factor that up scatters seed photons at 
the break to the desired observed energy. Eor the cases of inter¬ 
est, 6 <C 1 /cb, and the term in the radical dominates as shown. 
In the limit that e <C 76 , the remaining integrals may then be 
approximated by 
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where we have assumed that (5 Kljl, and 
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where in the final expression the factors of m7 have been rein¬ 
troduced. The typical flux at energy e is then 
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